Ultrafast  photovoltaic  response  in  ferroelectric  nanolayers 
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Abstract 


We  show  that  light  drives  large-amplitude  structural  changes  in  thin  films  of  the 
prototypical  ferroelectric  PbTiCfi  via  direct  coupling  to  its  intrinsic  photovoltaic  response.  Using 
time-resolved  x-ray  scattering  to  visualize  atomic  displacements  on  femtosecond  timescales, 
photoinduced  changes  in  the  unit-cell  tetragonality  are  observed.  These  are  driven  by  the  motion 
of  photogenerated  free  charges  within  the  ferroelectric  and  can  be  simply  explained  by  a  model 
including  both  shift  and  screening  currents,  associated  with  the  displacement  of  electrons  first 
antiparallel  to  and  then  parallel  to  the  ferroelectric  polarization  direction. 
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Text 


Light  couples  to  atomic-scale  degrees  of  freedom  in  complex  materials,  offering  new 
avenues  for  engineering  functionality  in  nanoelectronic  devices  and  enhancing  material 
properties.  The  key  functionality  of  ferroelectrics — a  thermodynamically  stable,  switchable 
polarization  that  persists  even  in  the  absence  of  an  applied  electric  field — itself  destabilizes  the 
ferroelectric  phase  in  thin  films,  due  to  the  existence  of  an  internal  depolarization  field  [1-7], 
associated  with  the  surface  charge  density  a  =  P  •  n  (where  P  is  the  electric  polarization  and  n  is 
the  unit  normal  vector),  that  points  opposite  to  the  direction  of  the  polarization.  This  field  can  be 
screened  by  surface  adsorbates  or  free  charges  [1,5],  as  in  the  simple  case  of  ferroelectric 
capacitors  with  metallic  electrodes,  and  can  lead  to  the  formation  of  stripe  domains:  periodic 
nanometer-scale  domains  of  alternating  polarization  that  minimize  the  free  energy  of  the 
system  [3,4,8].  Intensive  research  has  been  aimed  at  bypassing  the  intrinsic  size  limits  imposed 
by  the  depolarization  field  on  the  ferroelectric  phase,  including  epitaxial  strain 
engineering  [2,9,10]  and  the  use  of  electrodes  to  screen  the  depolarization  field  [6,1 1]. 

At  the  same  time,  ongoing  work  over  several  decades  [12-21]  has  sought  to  elucidate 
photovoltaic  effects  in  ferroelectrics  and  multiferroics  and  the  role  of  the  depolarization  field  in 
these  phenomena  [15,19].  These  are  associated  with  the  splitting  of  electron-hole  pairs  within  a 
noncentrosymmetric  crystal  and  lead  to  the  generation  of  open-circuit  voltages  much  larger  than 
the  bandgap.  A  number  of  mechanisms  have  been  implicated  in  these  effects,  including  the 
presence  of  the  internal  depolarization  field  in  thin  films  [15],  internal  fields  at  domain 
walls  [18],  and  second-order  nonlinear  optical  responses  in  the  field  of  the  incident  light  in  the 
medium  [21,22],  Such  nonlinear  responses  include  (non-resonant)  optical  rectification  [21];  shift 
currents  (also  referred  to  as  the  “linear  photovoltaic  effect”  [22]),  associated  with  an  anisotropic 
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displacement  of  electrons  and  holes;  and  injection  currents  [22],  Little  is  known  about  the 
coupling  of  these  photogenerated  currents  to  the  structural  response  of  ferroelectrics.  In  addition, 
the  possibility  of  using  light  to  drive  structural  changes  has  been  of  continuing  interest  [23-27]. 
Previous  time-resolved  x-ray  studies  have  focused  on  ferroelectric  heterostructures,  in  which 
photons  are  absorbed  within  metallic  layers  in  order  to  launch  acoustic  strain  responses  within 
the  ferroelectric  [27].  Here,  using  femtosecond  time-resolved  x-ray  scattering,  we  report  changes 
in  tetragonality  lasting  several  nanoseconds  in  the  prototypical  ferroelectric  perovskite  PbTiCfi 
(PTO),  following  photoexcitation  by  ultrafast  400  nm  light  pulses.  This  result  can  be  understood 
by  taking  into  account  the  coupling  of  the  photogenerated,  time-dependent  internal  fields  to  the 
strain  within  the  film. 

We  performed  femtosecond  time-resolved  x-ray  diffraction  measurements  of  a  20  nm 
PTO  film  on  SrTi03  (STO)  and  a  100  nm  PTO  film  on  DyScCfi  (DSO),  without  intermediary 
electrodes,  in  air,  under  photoexcitation  by  ultrafast  400  nm  laser  pulses.  Due  to  the  differing  in¬ 
plane  strains  provided  by  STO  (001)  and  the  pseudocubic  DSO  (110)  face,  Tc  of  the  two  films  is 
680°  C  and  470°  C,  respectively  [9].  We  studied  PTO  on  DSO  in  order  to  access  states  above  Tc 
while  avoiding  the  decomposition  of  the  film  at  high  temperature  due  to  loss  of  volatile  PbO  [8], 
Femtosecond  400  nm  laser  pulses  were  generated  via  second-hannonic  generation  from  800  nm, 
40  fs  full  width  at  half  max  (FWHM)  Tiisapphire  laser  pulses  and  incident  upon  the  samples  with 
5-polarization  at  a  60°  angle  from  the  nonnal,  with  a  maximum  incident  fluence  of  5  mJ/cm2 
(~  100  GW/cm2).  X-ray  pulses  with  60  fs  FWHM  [28]  and  X  =  1.39  A  were  scattered  from  the 
films  after  being  monochromatized  in  a  Si  (1 1 1)  double-crystal  monochromator  at  the  Linac 
Coherent  Light  Source,  giving  ~1010  photons/pulse  incident  on  the  sample  at  a  60  Hz  repetition 
rate.  We  observed  no  evidence  of  x-ray-induced  sample  modification.  The  time  resolution  of  this 
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measurement  is  limited  by  timing  jitter  between  laser  and  x-ray  pulses,  estimated  to  be 
100  fs  RMS  [29].  The  PTO  (003)  Bragg  reflection  was  collected  with  either  a  point  detector  or 
an  area  detector  as  a  function  of  laser  delay  and  x-ray  incidence  angle  (0),  probing  the  out-of¬ 
plane  lattice  parameter,  c,  inside  the  thin  film.  Similar  experiments  were  perfonned  at  the 
Advanced  Photon  Source,  but  temporal  resolution  was  limited  to  100  ps  by  the  synchrotron  x-ray 
pulse  duration.  Thin  films  of  PTO  were  deposited  on  STO  (001)  and  DSO  (110)  substrates  at 
930-990  K  via  metal-organic  chemical  vapor  deposition,  following  procedures  described 
elsewhere  [30],  Prior  to  growth,  Ti02-tenninated  surfaces  of  STO  were  created  via  a  standard  HF 
etch  [31].  DSO  surfaces  were  prepared  with  ScCF  termination  by  annealing  in  O2,  followed  by  a 
NaOH  etch  [32], 

When  grown  on  STO  (001),  PTO  forms  out-of-plane,  c-axis-oriented  domains  below  the 
Curie  temperature,  Tc,  due  to  compressive  strain  from  lattice  mismatch  to  the  STO  [2,3,9]. 

Under  zero-field  heating  conditions,  the  c-axis  of  PTO  contracts  as  Tc  is  approached, 
corresponding  to  a  negative  thermal  expansion  coefficient.  This  contraction  indicates  a  decrease 
in  the  tetragonality  and  the  dipole  moment  of  the  unit  cell,  and  causes  out-of-plane  diffraction 
peaks  to  shift  to  higher  scattering  angles.  In  contrast  to  the  expected  thermally  induced 
responses,  we  observe  significantly  more  complex  dynamics.  Time-dependent  rocking  curves  of 
photoexcited,  monodomain  PTO  on  STO,  with  polarization  pointing  out  of  the  film  [33],  taken  at 
room  temperature  [Fig.  1(a)  and  Fig.  1(b)]  indicate  a  symmetric  shift  of  the  Bragg  peak  to  higher 
Q  [where  Q  =  (47t  sin  0)/A]  within  5  ps.  This  is  consistent  with  a  uniform  contraction  of  the  unit 
cell  in  the  out-of-plane  direction  and  occurs  on  an  acoustically  limited  timescale,  determined  by 
the  film  thickness  over  the  sound  velocity  (y,  ~  4,000  m/s)  [34].  Subsequently,  we  observe  a 
long-lived,  symmetric  shift  of  the  diffraction  peak  to  lower  Q  occurring  on  a  10  ps  timescale,  in 
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contrast  to  what  would  be  expected  from  simple  heating  of  a  ferroelectric.  Following  the  low-Q 
shift,  the  nanolayer  relaxes  to  its  initial  state  before  excitation  on  single-nanosecond  timescales 
[Fig.  1(c)].  We  note  that  the  in-plane  lattice  constants  are  expected  to  evolve  on  acoustically 
limited  timescales  detennined  by  the  1-inm  laser  spot  size,  and  therefore  do  not  change  on  the 
picosecond  timescales  here. 

In  order  to  further  elucidate  the  time-dependent  changes  in  the  ferroelectric  unit  cell,  we 
measured  the  response  of  the  system  at  535°  C,  where  PTO  on  STO  enters  a  stripe  domain 
phase  [3],  characterized  by  neighboring  regions  of  opposite  polarization,  with  a  period  A  of  order 
10  mn.  This  period  is  detennined  by  a  trade-off  between  the  energy  of  the  uncompensated 
depolarization  field  and  the  domain  walls  [8],  At  higher  temperatures,  the  volatile  adsorbate 
species  that  stabilize  the  monodomain  phase  [1,5]  are  removed.  The  periodicity  of  the  stripe 
domains  gives  rise  to  in-plane  satellites,  displaced  by  AQr=27t/A,  which  appear  as  a  ring  of 
diffuse  scattering  around  the  Bragg  peak.  We  measured  the  time  evolution  of  this  diffuse 
scattering  under  photoexcitation  (Fig.  2).  After  correcting  the  orientation  of  the  imaging  plane  to 
make  it  nonnal  to  the  Z-axis  in  reciprocal  space  [Fig.  2(a)],  the  distance  in  reciprocal  space  from 
the  Bragg  peak  to  the  center  of  the  diffuse  scattering  was  defined  to  be  Qr  =  (Qh2  +  Qk2)12-  Ql 
was  defined  as  the  out-of-plane  displacement  of  the  diffuse  scattering.  The  center  of  mass  of  the 
diffracted  intensity  was  then  determined  as  a  function  of  laser  delay  [Fig.  2(b)].  We  observe  a 
transient  decrease  in  the  out-of-plane  lattice  constant  (measured  by  Ql)  and  the  stripe  domain 
period  (measured  by  Qr),  followed  by  a  long-lived  increase  on  a  timescale  matching  the 
monodomain  room  temperature  response.  The  strong  match  between  the  timescales  in  both 
phases  [Fig.  2(b)]  indicates  a  mechanism  independent  of  the  intrinsic  domain  structure  and  rules 
out  effects  associated  with  photothermal  desorption. 
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Based  on  the  measurements  described  above,  we  propose  the  following  two-step 
mechanism  for  the  observed  dynamics,  driven  by  the  displacement  of  electrons  first  against  and 
then  along  the  polarization  direction  (pointing  out  of  the  film)  (Fig.  3).  In  the  first  step,  the  initial 
decrease  in  the  ferroelectric  tetragonality  during  the  first  5  ps  follows  from  a  shift  current, 
associated  with  a  photoinduced,  coherent  shift  of  electrons  away  from  the  free  surface  [35,36]. 
This  current,  directed  along  the  polarization  direction,  increases  the  charge  density  at  the 
interfaces  [37],  giving  rise  to  an  increase  in  the  depolarization  field  and  an  impulsive  stress 
which  leads  to  contraction  of  the  lattice.  We  estimate  this  contribution  by  calculating  the  second- 
order  nonlinear  current  in  the  applied  light  field,  obtained  by  a  first-principles  calculation  of  the 
third-rank  shift  current  tensor,  in  order  to  detennine  the  photoinduced  current  density  while  the 
sample  is  illuminated  [38].  For  PTO  (point  group  4 mm),  the  relevant  tensor  components  have 
subscripts  31,  33,  and  15;  the  31  and  33  components  represent  the  current  parallel  to  the 
ferroelectric  polarization;  and  their  values  depend  on  the  light  frequency  and 
polarization  [21,38].  Our  calculation  of  the  current  direction  is  consistent  with  previous 
observations  [13].  Quantitatively,  we  estimate  surface  charge  densities  of  5  flC/cm2  for  the 
fluences,  light  polarization  and  excitation  wavelength  used  here.  In  the  second  step,  after  the 
light  is  turned  off,  the  transient  current  reverses  direction  as  the  generated  electrons  and  holes 
move  in  and  screen  the  depolarization  field,  leading  to  an  increase  in  the  c-axis  lattice  parameter. 
This  model  is  supported  by  several  key  observations: 

(1)  As  noted  before,  concomitant  with  the  changes  in  tetragonality,  we  observe  dynamical 
changes  in  the  stripe  domain  period,  A,  consistent  with  modification  of  the  surface  charge  density 
and  the  depolarization  field.  Increases  (decreases)  in  A  reflect  more  complete  (less  complete) 
screening  of  the  depolarization  field  by,  for  example,  compensating  ions,  surface  defects,  or  free 
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charges  [1,4, 5, 7].  These  changes  are  analogous  to  the  compensation  effects  observed  in 
ferroelectrics  with  electrode  structures  [7], 

(2)  Corresponding  measurements  of  PTO  on  DSO  in  the  paraelectric  phase  [Fig.  4(a)], 
where  the  depolarization  field  vanishes  due  to  the  centrosymmetric  unit  cell,  show  negligible 
shifts  in  0  (<0.004°)  during  the  first  20  ps.  The  centrosymmetry  of  the  paraelectric  phase  also 

precludes  second-order  nonlinear  optical  responses  such  as  shift  currents,  in  agreement  with  our 
observations. 

(3)  In  Fig.  4(b),  the  maximum  observed  c-axis  lattice  parameter  for  PTO  on  STO  as  a 
function  of  temperature  is  compared  to  the  c-axis  parameter  without  a  depolarization  field 
calculated  from  a  Landau-Ginzburg-Devonshire  (LGD)  model,  which  predicts  the  equilibrium 
relationships  between  the  polarization,  c-axis  parameter,  internal  field,  and  epitaxial  strain  seen 
by  the  PTO  [9,39].  The  fully  compensated  temperature-dependent  polarization  and  resulting  c- 
axis  parameter  were  calculated  by  minimizing  the  bulk  LGD  free  energy  density  under  zero 
internal  electric  field  (i.e.,  zero  depolarization  field)  conditions  [9],  describing  the  case  when  the 
free  carriers  fully  compensate  the  interfacial  charge.  This  model  considers  the  effect  of  epitaxial 
strain  on  the  thermodynamic  stability  of  ferroelectric  PTO  phases  and  accounts  for  the  known 
temperature  dependence  of  the  STO  lattice  parameter  [40].  Thermodynamic  constants  relevant  to 
the  ferroelectric  properties  of  PTO  were  taken  from  literature  values  [41],  These  numbers  are  in 
good  agreement  and  are  consistent  with  nearly  complete  screening  of  the  internal  field. 

(4)  As  shown  in  Fig.  4(c),  the  induced  strain  at  t  =  100  ps  saturates  with  increasing 
fluence.  This  observation  is  consistent  with  a  free-carrier-induced  screening  of  the  internal  field 
which  saturates  when  the  field  is  fully  canceled  [13,42,43]  and  is  inconsistent  with  models 
associated  with  pyroelectric  or  other  temperature-induced  effects,  which  would  result  in  c-axis 
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lattice  parameter  changes  in  the  opposite  direction  at  long  times  and  are  not  expected  to  saturate. 
We  measure  an  absorption  constant  apio  =  1  x  103  cm"1  at  X  =  400  mn  by  ellipsometry  (in 
agreement  with  our  first-principles  calculations  of  the  absorption  coefficient  [38]),  from  which 
we  estimate  a  carrier  density  N  ~  1020/cm3  and  a  temperature  jump  of  approximately  50  K  [44] 
for  the  fluences  used  here. 

The  transient  polarization  change,  A Peh,  associated  with  an  electron-hole  displacement  Az, 
is  A Peh  =  NqAz  ~  10  p,C/cm2  for  Az  comparable  to  the  nanolayer  thickness,  consistent  with  the 
shift-current-induced  modulation  calculated  above  [45].  We  scale  the  maximum  observed  Ac  at 
t  >  10  ps  using  the  equilibrium  electrostrictive  response  Ac  ~  \P\2  [39]  to  estimate  the  magnitude 
of  the  ferroelectric  polarization  increase,  and  obtain  roughly  5  ptC/cm2,  a  significant  fraction  of 
the  initial  polarization  (Pq  ~  60  p. C/cm 2  [2]).  Carrier  motion  in  the  internal  field  of  the  film  is 
expected  to  develop  on  a  timescale  of  roughly  A z/(jiE)  ~  5  ps,  using  typical  values  for  the  carrier 
mobility  (/t  ~  1  cm2  V"1  s"1)  [46]  and  the  internal  field  ( E  ~  500  kV/cm)  [39]  in  PTO.  It  is 
therefore  delayed  with  respect  to  the  shift  current  response,  in  agreement  with  our  observations. 
At  all  temperatures,  only  small  changes  in  the  STO  or  DSO  diffraction  peaks  are  observed  (not 
shown).  Measurements  using  well-above-bandgap  266  mn  pulses  (not  shown)  show  effects  of 
similar  magnitude  at  significantly  lower  fluences,  of  order  10  pJ/cm^  consistent  with  the 
increased  optical  cross  section  at  266  mn.  (We  measure  an  absorption  constant 
ocpio  =  9  x  105  cm"1  at  X  =  266  mn.)  This  observation  precludes  optical  rectification,  which  is 
expected  to  depend  primarily  on  incident  intensity,  as  a  possible  mechanism  for  these  effects, 
consistent  with  previous  work  [22]. 
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To  further  interpret  Fig.  1,  we  model  the  strain  evolution  within  the  PTO  as  a  response  to 
an  applied  stress  by  the  photogenerated  currents  [47].  The  acoustic  wave  equation  is  solved  for  a 
thin  film-substrate  system  for  a  time-dependent  pressure  described  by  the  following  equation: 

a(t)  =  H(t)  { a  sc  +  oDP  [1  -  exp(-//r)]  ]  (Eq.  1) 

Here,  H(t)  is  the  Heaviside  step  function  and  c  sc  represents  an  instantaneous  stress  associated 
with  the  shift-current-induced  increase  in  the  depolarization  field.  The  subsequent  tetragonality 
increase,  ascribed  to  screening  of  the  internal  field  by  free  carriers,  leads  to  a  slowly  increasing 
stress,  <5 dp,  in  the  opposite  direction,  with  a  time  constant  r  determined  by  the  carrier  mobility. 
These  stresses  develop  from  the  piezoelectric  coupling  between  the  lattice  and  the  internal  field. 
Transient  PTO  (003)  rocking  curves  resulting  from  the  stress  described  by  Eq.  1  are  calculated 
using  dynamical  diffraction  theory  [48]  and  used  to  construct  a  predicted  time  scan  like  that  of 
Fig.  1.  As  can  be  seen  by  comparison  to  the  measured  data  (Fig.  3),  this  simple  acoustic  model 
provides  an  excellent  description  of  the  experimental  results  during  the  first  tens  of  picoseconds. 
The  best  quantitative  fitting  is  obtained  for  Odp/Osc  =  -1.6  and  r=  4  ps,  in  good  agreement  with 
the  estimates  given  above  for  the  screening  time. 

The  strong  coupling  between  light,  photocurrents,  and  atomic-scale  degrees  of  freedom 
elucidated  here  determine  the  first  steps  in  the  processes  that  underlie  the  intrinsic  photovoltaic 
response  of  ferroelectrics,  and  result  from  a  current-driven  modulation  of  the  internal  field  within 
the  thin  film.  This  work  provides  new  avenues  for  enhancing  the  functionality  of  ferroelectric 
devices  on  ultrafast  timescales  by  a  light-mediated  coupling  to  their  polarization.  Further 
engineering  of  artificial  ferroelectric  structures  and  extensions  of  this  work  to  light  frequencies 
spanning  the  range  from  visible  to  terahertz  may  soon  provide  a  technologically  viable  and  facile 
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pathway  to  ultrafast  all-optical  switching,  addressing  a  key  obstruction  towards  the  goal  of 
memory  devices  both  written  to  and  read  out  by  light. 
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Figures 

FIG.  1.  Evolution  of  the  tetragonality  of  PTO  on  STO  in  the  monodomain  phase. 

(a)  Time  scans  on  the  low-  and  high-angle  sides  of  the  PTO  (003)  rocking  curve  at  room 
temperature  (incident  fluence  =  5  mJ/cm2).  We  observe  an  acoustically  limited  high-Q  shift, 


followed  by  a  long-lived  low-Q  shift,  (b)  Rocking  curve  scans  at  various  time  delays  denoted  in 
(a),  (c)  Long-range  time  scan  on  the  low-angle  side  of  PTO  (004)  at  350°  C  (incident  fluence  = 
5  mJ/cm2).  The  observed  low-Q  shift  decays  on  single-nanosecond  timescales. 
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FIG.  2.  Stripe  domain  phase  response  at  535°  C.  (a)  Area  detector  images  of  PTO  (003) 
at  535°  C  for  several  time  delays  (incident  fluence  =  5  mJ/cm2),  corrected  such  that  the  angle  of 
the  detector  plane  is  normal  to  the  /.-axis  in  reciprocal  space,  (b)  Time  evolution  of  Ql  and  Qr, 
extracted  from  area  detector  images  of  PTO  (003).  Ql  is  the  position  of  the  diffuse  scattering 
ring  along  the  /.-axis  and  Qr  is  the  radius  of  the  diffuse  ring. 
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FIG.  3.  Schematic  of  unit  cell  response.  At  t  =  5  ps,  the  c-axis  lattice  parameter  reaches  a 
minimum  due  to  the  piezoelectric  response  to  the  increase  in  the  depolarization  field,  Ed,  induced 
by  the  shift  current.  At  t  >  10  ps,  carriers  have  moved  to  screen  the  depolarization  field,  driving  a 
long-lived  increase  in  tetragonality.  Time  scans  from  Fig.  1  are  fit  by  modeling  the  out-of-plane 
strain  of  PTO  as  a  response  to  the  time-dependent  stress  profile  given  in  the  inset  (see  Eq.  1). 
Vertical  lines  are  guides  to  the  eye. 
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FIG.  4.  Role  of  depolarization  field  in  photovoltaic  response,  (a)  PTO  (003)  on  DSO  in 
the  paraelectric  phase  at  479°  C,  slightly  above  Tc  (incident  fluence  =  5  mJ/cm2).  Only  small  0 
shifts  are  seen  due  to  the  centrosymmetric,  paraelectric  unit  cell,  (b)  LGD  theory  calculation, 
comparing  the  observed  c-axis  lattice  parameter  of  PTO  on  STO  at  t  <  to  and  t  =  100  ps  (incident 
fluence  =  4  mJ/cm2)  to  the  predicted  zero-field  c-axis  parameter,  (c)  Strain  in  the  c-axis 
parameter  of  monodomain  PTO  on  STO  at  550°  C  as  a  function  of  fluence,  measured  at  t  =  100 
ps.  Saturation  behavior  in  the  maximum  drivable  strain  is  seen. 
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